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ABSTRACT 

I present a general discussion of the evolution and model-dependence of both the 
linear and nonlinear power spectrum of density fluctuations. The features of the 
linear power spectrum in cosmological models with cold dark matter (CDM) and 
cold+hot dark matter (C+HDM) are reviewed, and useful analytical approximations 
are summarized. Cosmological numerical simulation results are then used to illustrate 
the process of nonlinear gravitational clustering and to compute the nonlinear power 
spectrum. A new analytical approximation relating the linear and nonlinear power 
spectrum in C+HDM models is presented. 

1. Introduction 

The power spectrum P(k) of matter fluctuations is a fundamental quantity in cosmology. 
It provides the most basic statistical measure of gravitational clustering, and an accurate 
determination of the power spectrum is among the most important goals of every major galaxy 
survey. Since each cosmological model has its distinct power spectrum, one can hope to obtain 
crucial information such as the matter content of the Universe and the nature of the primordial 
fluctuations by comparing the observed power spectrum with theoretical predictions for various 
plausible models. 

Measurements of the power spectrum require extensive galaxy surveys covering a large 
volume of space. It is only recently that such surveys have been completed and are starting to 
reveal intriguing results. For instance, much interest has been generated by the compilation 
of the power spectrum determined from the APM and IRAS galaxy catalogs and the effort to 
reconstruct the underlying linear mass power spectrum (Peacock & Dodds 1994; Peacock 1997). 
The recently completed Las Campanas survey of over 23,000 galaxies has yielded a power spectrum 
incompatible with the standard cold dark matter (CDM) model but still consistent with a class of 
open CDM, cold+hot dark matter (C+HDM), and tilted CDM models (Lin et al. 1996). Larger 
ongoing projects such as the CfA Century, the 2dF, and the Sloan surveys promise to provide 
measurements with reduced error bars and data points on larger scales for better constraints. 

On the theoretical side, calculations of the power spectrum falls into two regimes: the 
linear and the nonlinear, which are characterized by the amplitude of the density fluctuations 
5 = 5p/p <C 1 and > 1, respectively. Since the observed power spectrum spans both linear and 
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nonlinear regimes, theoretical calculations in both regimes must be performed before all data points 
can be fully utilized as constraints. In the linear regime, the computation of the power spectrum 
has become standard practice. It is obtained by integrating the coupled, linearized Einstein 
and Boltzmann equations that describe the evolution of the metric and density perturbations 
throughout the history of the Universe. Several groups have written numerical codes for such 
calculations, and a publicly-available version is described in Ma & Bertschinger (1995) and placed 



at http://arcturus.mit.edu . This code computes, in both synchronous and conformal Newtonian 
gauges, the evolution of the phase-space distributions of photons, baryons, cold dark matter, and 
both massless and massive neutrinos. It is therefore applicable to most CDM, C+HDM, and 
CDM with a cosmological constant (LCDM) models. To gain deeper insight and for convenience, 
analytical approximations for the linear power spectra in these models have also been published. 
Details about the linear power spectrum are discussed in Section 2. 

Gravitational clustering, however, eventually becomes a nonlinear process. The existence 
of galaxies and clusters is a manifestation of the nonlinear nature of gravity. Determining the 
nonlinear power spectrum P n \ is therefore an important task. Not surprisingly, P n \ is more difficult 
to obtain than the linear power spectrum. Higher-order perturbation theories help to extend 
the range of validity to the quasi-linear regime, but the fully nonlinear power spectrum can be 
calculated only from numerical simulations. Thus far, theoretical predictions of P n \ have been 
carried out for only a few models at limited epochs because it is difficult and time-consuming 
to perform numerical simulations with sufficient dynamic range to allow calculations of P n \ over 
a wide range of scales. Section 3 summarizes some recent work and presents new analytic and 
numerical results for nonlinear clustering in C+HDM models. 



2. Linear Power Spectrum 

The power spectrum P(k) of matter fluctuations is related to the density field 6 in fc-space by 

(S(k)6*(k')) = P(k)5 D (k - k') , (1) 

where 5d is the Dirac-delta function. The Fourier transform of the power spectrum is the two-point 
correlation function £(r) = (5(x{)5(x2)) , where r = \x\ — x%\. For a Gaussian density field (as 
predicted in most inflationary theories), its statistical property is entirely determined by P(k). 

The shape of the linear power spectrum depends on the cosmological parameters assumed 
in a given theory of structure formation. It is governed by the linear perturbation theory of 
gravitational clustering, and can be computed by integrating the coupled, linearized Einstein, 
Boltzmann, and fluid equations for the metric and density perturbations. Some important 
parameters that affect the shape of the power spectrum are: (1) Primordial spectral index n: 
P(k) oc k n . An example is the Harrison-Zeldovich (Harrison 1970; Zeldovich 1972) spectrum 
which takes n = 1. Most versions of inflationary models also predict nearly n = 1 power spectra. 
(2) Matter-radiation equality time: The equality time i eq is defined to be the epoch in the thermal 
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history of the Universe when the energy density in matter equals that in radiation. The equality 
redshift scales as 1 + z eq oc tl m h 2 , where fi m is the density parameter in matter and h is the Hubble 
constant in units of 100 km/s/Mpc. This parameter controls the location of the peak of P(k) in 
CDM-type models since the density perturbations with wavenumbers k > k eq enter the horizon in 
the radiation-dominated era and cannot grow appreciably, whereas perturbations with k < k cq 
enter the horizon in the matter-dominated era and can grow unimpeded. (3) Nature of dark 
matter. The pure CDM model, for example, exhibits a characteristic P(k) oc A;~ 3 slope at high 
k, while the power spectrum for the pure hot dark matter (HDM) model is cut off exponentially 
below the free-streaming scale of the massive neutrinos due to the phase mixing in the neutrino 
phase-space distribution. The hybrid C+HDM models, which are parameterized by the neutrino 
mass density Q. v and CDM mass density O c = 1 — Vt p — Vt^ (where fi& is the mass density in 
baryons), exhibit intermediate behavior. The suppression in the power in both the cold and hot 
components due to neutrino free-streaming generally increases with increasing Q. v . 

The solid curves in Figure 1 show the linear power spectra for the pure CDM and three 
C+HDM models with £l u = 0.1,0.2, and 0.3. The corresponding neutrino masses in the four 
models are 0, 2.3, 4.6, and 7 eV, respectively. (Only one of the three species of neutrinos is 
assumed massive.) The models have a total matter density of Q, m = 1 and h = 0.5, and all are 
normalized to the COBE rms quadrupole Q rms -ps = 18 fiK (Gorski et al. 1996), as evidenced by 
the convergence of the curves at small k. The differences at large k reflect the distinct properties 
of cold and hot dark matter: the larger Q, v is, the weaker is the gravitational clustering on 
small length scales. For comparison, a COBE-normalized low-density model with Q m = 0.3, a 
cosmological constant fl\ = 0.7, and h = 0.75 is also shown (dashed curve). 

Analytical approximations with less than 10% error are very convenient as an input for many 
calculations. The pure CDM model is well-approximated by 



where A is the overall normalization factor, a is the expansion factor, n is the primordial spectral 
index, k is the wavenumber in units of Mpc -1 , and a± = 2.34, 0:2 = 3.89, «3 = 16.1, = 5.46, 
and 05 = 6.71 (Bardeen et al. 1986). The normalization A can be determined from the 4- year 
data from the COBE satellite, and A = 2420 /i -4 Mpc 4 for n = 1 flat models (Bunn & White 
1997). The shape parameter T characterizes the dependence on cosmological parameters, and 
a good approximation is found to be given by F = tt m h/ exp[f2t,(l + l/O m )] (Efstathiou et 
al. 1992; Sugiyama 1995). The error in the approximation for the standard CDM model with 
Ob = 0.05, for example, is < 10%. A higher-accuracy fit for this much-studied model can be 
achieved by setting Ob = in equation @ (i.e. setting T = Q m h) and modifying the coefficients 
to a\ = 2.205, «2 = 4.05, «3 = 18.3, 04 = 8.725, and 05 = 8.0. The fractional error relative to the 
direct numerical result is then reduced to smaller than 1% for k < 40 /i Mpc -1 (Ma 1996). 
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The linear power spectra for the C+HDM models require additional treatment since the effect 
of neutrino free-streaming on the shape of the power spectrum is both time- and scale-dependent. 
It is found (Ma 1996) that by introducing a second shape parameter 

T u = a l / 2 tt u h 2 (3) 

to characterize the neutrino free-streaming distance, accurate approximations to the linear power 
spectra in C+HDM models can be obtained. The cold and hot components clearly have different 
power spectra due to their different thermal properties. For the CDM component in C+HDM 
models, a good approximation is given by 

(. Ql.05 
I _|_ }) l x b4 / 2 + 6 2 x bA \ k 
bl > x = ?r > xo=x(a = l), (4) 
l + b 3 x 4 J U 

where P c (k, a,£l u = 0) for the pure CDM model is given by equation (Q), and 
h = 0.01647, b 2 = 2.803 x lO" 5 ,^ = 10.90, and b 4 = 3.259 for k in units of Mpc" 1 . The 
functional form for the ratio P c (k, a, £l u )/P c (k, a, 0) is chosen to have the asymptotic behavior 
oc a 2 ^°°~ l \ which can be derived analytically, and /oo = (5-^/1 — 24^/25 — l)/4 is the asymptotic 
growth rate for k — > oo, where 1 — /oo oc f}J/ 05 is a good approximation. For the HDM component 
in C+HDM models, an accurate approximation is given by 

P u (k,a,n u ) = P c (k,a,n u ) [ — 6 ■ _ -) , x = -\ , (5) 

1 + C2 X ' + C3 X + C4 X A ' Z + C5 X z I V v h 



where a = 0.0015, c 2 = -0.1207, c 3 = 0.1015, c 4 = -0.01618, and c 5 = 0.001711 for k in units of 
Mpc" 1 . 

It is often useful to have an analytic approximation for the density-weighted power spectrum 
p(k) = {n u Pu /2 + (i - n u )Pc /2 } 2 that measures the total gravitational fluctuations contributed 
by the two components. Here the CDM and baryons have been assumed to have the same power 
(i.e., P c = Pj,), which is a good approximation for the range of redshifts and studied in this 
paper. The functional form used for the CDM spectrum P c in equation (Q) works well here, and a 
good approximation for the density-weighted power spectrum in C+HDM models is given by 

P(k,a,Sl v ) , _ . / l+ji x d */ 2 + d 2 x d * \ nl '° 5 k , n ((Cs 

- g(x, tt u ) = — -r , x = — , x = x[a = 1) , (6) 



p c (k,a,si v = o) K ' v i + rfs^ 4 J ' r 

where P c (k,a, Q u = 0) again is given by equation (2), and the coefficients are 
d x = 0.004321, d 2 = 2.217 x 10~ 6 , d 3 = 11.63, and d 4 = 3.317 for k in units of Mpc" 1 . 

It should be noted that the C+HDM power spectra do not obey the simple evolution 
P(k) oc a 2 in equation (2) for the flat pure CDM model. The free-streaming of massive neutrinos 
slows down with time, allowing the neutrinos to cluster gravitationally and the neutrino density 
perturbations to grow on increasingly smaller scales. As a result, the growth of the C+HDM power 
spectra is both scale- and time-dependent. This effect is taken care of by the time-dependent 
parameter T u in equation (3), which is built in in equations (4)- (6) via the scaled variable x. 
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3. Nonlinear Gravitational Clustering 

3.1. Numerical Simulations 

The linear theory of gravitational clustering is an elegant and powerful theory. It describes 
accurately the growth of density perturbations from the early Universe until a redshift of ~ 100, 
and the solution to the theory involves straightforward time-integration of a set (albeit a large set) 
of coupled ordinary differential equations. Gravitational clustering, however, eventually becomes 
a nonlinear process, and the study of the fully nonlinear process ultimately relies on numerical 
simulations. 

Cosmological simulations can be roughly divided into two categories: (1) iV-body 
simulations that deal only with dissipationless gravitational interactions of dark matter; and (2) 
hydrodynamical simulations that model gaseous dissipation via cooling and heating processes 
in addition to the gravitational interactions among the dark matter and gas. Depending on 
the models, the simulations are started at a redshift between 20 and 100 when the rms density 
fluctuations are well below unity, and the initial conditions of the simulations are generated from 
the linear power spectrum at that epoch. 

As an illustration of gravitational clustering in collisionless iV-body simulations, Figure 2 
shows the projection of the smoothed matter distribution at four redshifts, z = 5,3, 1.5, and 0, 
from an iV-body simulation of the VL V = 0.1 C+HDM model. The model assumes £l c = 0.85, 
f&b = 0.05, and h = 0.5, and is normalized to the COBE quadrupole Qrms-PS = 18 fiK (Gorski 
et al. 1996). The simulation box is 100 Mpc comoving on a side, and the force softening is 50 
kpc. A total of 128 3 cold and 128 3 hot particles are used. Each panel in Figure 2 is 100x100 
Mpc comoving, and shows the projection along one axis of the entire simulation box. Periodic 
boundary conditions are adopted in the simulation, so for example, the dense clumps near the 
central left and right edges of the box in the z = panel (lower right) are part of the same cluster 
of dark matter halos. The darkest several halos at z = all have masses above 1O 14 M . 

The power spectrum of the density fluctuations offers the lowest-order measurement of the 
gravitational clustering exhibited in Figure 2. In Figure 3, we plot the corresponding nonlinear 
power spectra computed from the particle spatial distributions at z = 3, 1.5, and shown in 
Figure 2. The linear power spectra given by equation (6) are also shown for comparison. The 
hierarchical, or "bottom-up", nature of gravitational collapse in these models is evident: the 
high-/? modes (i.e. small length scales) have become strongly nonlinear, while the low-A: modes are 
still following the linear power spectrum. The fact that the three lowest k modes are still linear at 
z = ensures that our choice of the simulation box size (100 Mpc) is large enough to include all 
waves that have gone nonlinear at the present epoch. It is also interesting to note that the point 
of departure from linearity moves to the left of the figure as z decreases, indicating that objects 
become nonlinear on increasingly larger length scales as the Universe evolves. 

It is also instructive to compare the nonlinear structures at the same epoch in different 
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cosmological models. Figure 4 shows the smoothed matter distribution at z = 1.5 from simulations 
of the four models in Figure 1: a LCDM, and C+HDM with VL V = 0.1,0.2, and 0.3 (clockwise 
starting from upper left). The LCDM model has fl m = 0.3, £l\ = 0.7, and h = 0.75, and 
the three C+HDM models have h = 0.5, f2b = 0.05, and f2 c + fi„ + Ob = 1- All models are 
COBE-normalized, and all simulations are performed with the same box size and force resolution 
as the VL V = 0.1 run described earlier. (The only exception is that the Q u = 0.2 and 0.3 simulations 
used 10 x 128 3 instead of 128 3 particles to represent the HDM component.) The initial conditions 
for all four simulations are generated with the same random phases; structures therefore appear 
at similar locations in all panels. The effect of £l v on structure formation, however, is striking: 
the higher Q, v is, the fewer collapsed objects there are at a given epoch. This trend simply 
reflects the decrease in the small-scale linear power shown in Figure 1, which results from neutrino 
free-streaming. Note that the model with Cl v = 0.3 (lower left panel), which corresponds to a 
neutrino mass of 7 eV, has very few structures even at z = 1.5. It is also interesting to note that 
the upper two panels of Figure 4 for the LCDM and fi„ = 0.1 C+HDM models look very similar. 
This is because the linear power spectra for the two models are in fact very similar at k > 0.1 h 
Mpc -1 , as shown in Figure 1. The corresponding linear and nonlinear power spectra at z = 1.5 
for the three C+HDM models in Figure 4 are shown in Figure 5. 

The power spectra shown thus far for the C+HDM models are the density-averaged 
P{k) = {n v Pu /2 + (1 - Q„)Pc /2 } 2 , where P v and P c are the individual power spectra for the 
hot and cold components. (It has been assumed P c = -Pbaryon, which is a good approximation for 
the redshift range of interest here.) The two components evolve distinctly due to their different 
thermal velocities, and the details of the shape and the growth of the linear power spectra are 
discussed in Ma (1996). Here, the nonlinear power spectra for the separate components are 
presented in Figure 6. For clarity, only two epochs, z = 3 and 0, are plotted. Notice how at z = 3, 
the HDM spectrum remains linear to almost k = 2/iMpc -1 but the CDM remains linear to only 
k ~ 0.5/iMpc" 1 , while at z = 0, both components become nonlinear at k ps 0.3 h Mpc -1 . This is 
largely due to the slowing down of the neutrinos, which makes it easier for the neutrino particles 
to fall into the CDM potential wells at later times, and therefore become nonlinear at smaller k. 
More precisely, the median thermal velocity of neutrinos is 

VRS 3(1 + z )Mk = 15(1 + z )H^X km / S) ( 7 ) 
m v c m v 

so the neutrinos in the 0„ = 0.1 model (m„ = 2.3 eV) have slowed down by a factor of 4 since 
z = 3 to about 60 km/s today. 



3.2. Analytical Approximations 

Just as it was useful to cast the linear power spectra in simple functional forms for a wide 
range of models (see equations (l)-(6)), analytical approximations to the nonlinear power spectrum 
can also provide valuable insight to the process of gravitational clustering. It is, however, more 
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difficult to compute the nonlinear P n \ because the numerical simulations that have sufficient 
dynamic range to allow calculations of P n \ over a wide range of scales are generally much more 
CPU-intensive than the integration of the linearized Boltzmann equations. The behavior of the 
power spectrum in the nonlinear regime therefore is less well-understood. 

An early attempt to relate linear and nonlinear quantities was focused on the spatially- 
averaged two-point correlation function £(r) in models with f2 m = 1 and a power-law power 
spectrum (Hamilton et al. 1991). It was found that if a given nonlinear scale r n \ is identified with 
its pre-collapsed linear scale r\ by rf = r^(l + £ n i), then there exists a simple, universal function 
relating the linear and nonlinear two-point correlation function: 

Ur„i) = FUn)] . (8) 

At that time, this transformation appeared to be magically independent of the spectral index n 
of the linear power spectrum assumed in the model, and the generality of this formula rendered 
the task of reconstructing the primordial spectrum from the observed nonlinear clustering of 
galaxies less complicated. However, further tests against numerical simulations showed that this 
formula worked well only for spectral index n > — 1, and it erred by factors up to 3 and 10 for 
n = —1.5 and —2, respectively (Jain et al. 1995). In particular, the nonlinear £ n i was found to 
rise more sharply with increasing £i for models with more negative n, so Jain et al. introduced an 
n-dependent formula to accommodate this feature. For the more realistic cosmological models such 
as the CDM where the spectral index is a function of scale, they proposed to use an effective index 
n e g, defined to be n e ff = din P/dln ko, where ko is the scale at which the rms mass fluctuation 
a(R = 1/ko) is unity. The index n c g therefore represents the slope of the power spectrum at a 
scale where nonlinear effects are becoming important. Later work extended the analysis to models 
with varying matter density Q m and cosmological constant Q\ (Peacock &: Dodds 1996), where 
instead of using an effective index for all scales, the local slope n{k) = din P/dln k was adopted. 

Here we examine the nonlinear mapping of the power spectrum in the C+HDM models, which 
has not been explored yet. Figure 7 shows how the nonlinear density variance A n i = ^nk^ P n i(k n \) 
diverges from the linear Ai = Airkf P\(k\), where the linear and nonlinear wavenumbers are related 
by kf = k^/(l + A n i). The squares are obtained from simulations of the £l v = 0.2 (left curve) and 
0.1 (right) C+HDM models. The curves generally obey the asymptotic condition A n i = Ai for 
Ai <C 1, and A n i oc Aj' 5 in the highly-nonlinear, stable clustering regime. The mapping, however, 
is clearly not independent of cosmology: the larger Q v is in a model, the faster A n i/Ai increases 
at Ai ~ 1. This trend reflects the n-dependence pointed out by Jain et al. (1995) and can be 
explained by the different shapes of Pi shown in Figure 1, where the models with higher 0, u have 
less power and hence more negative slope at high k. Since all models are normalized to COBE and 
have the same amplitude at low k, the critical ko (where cr(ko) = 1) increases for larger Q. v . The 
effective index n e g defined at ko is therefore more negative for larger Q v , resulting in the steeper 
rise of A n i in Figure 7. 

However, neither formula proposed by Jain et al. (1995) or Peacock & Dodds (1996) can be 
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extended to the C+HDM models. The dotted curves in Figure 7 illustrate the large discrepancies 
in the Peacock-Dodds fitting function when it is applied to the Cl u = 0.1 and 0.2 models. In 
retrospect, it is not surprising their formulas do not apply: Although the power spectrum at a 
given epoch is determined entirely from the spatial distribution of the particles, the evolution of 
the power spectrum depends on the particle velocities as well. A general formula for the linear to 
nonlinear transformation therefore must depend on both the shape and the growth rate of P(k). 
Since both formulas are designed for models without massive neutrinos and depend only on the 
shape of P(k), they cannot be applied to C+HDM models. 

Here we propose a new analytical approximation for the mapping of linear and nonlinear 
power spectrum in CDM as well as C+HDM models: f 



A nl (fc nl ) = Ai(fc 1 )G[A,(fci)] , (9) 

where 

_ 1 + Ql x 4 + a 2 x 8 /g 2 - 5 
G{X) - l + a 3 x* + a A x™/g*> ( } 

and the coefficients are a\ = 4756, a 2 = 384.6, = 3732, and = 24.20. The functional form of 
G is chosen to give the appropriate asymptotic behavior A n i — > Ai in the linear regime (x -C 1) 
and A n i oc Aj in the stable clustering regime (i > 1). The dependence of G on $7„ comes from 
the function g of equation (6), which gives the relative amplitudes of the power spectra in the 
C+HDM models and the pure CDM model. This function is analogous to the commonly-used 
growth factor g(a = 1, O m , A ) = f^m^m 7 - A + (1 + O m /2)(l + Oa/70)]" 1 for LCDM models 
(Lahav et al. 1991; Carroll, Press & Turner 1992). They differ, however, in that the growth factor 
is scale-independent in LCDM models but is a function of scale in C+HDM models since neutrinos 
only affect the growth below the free-streaming scale, as discussed in Section 2. 



4. Summary 

The power spectrum is a fundamental measure of gravitational clustering in cosmology. We 
have discussed the features and evolution of both the linear and nonlinear power spectra in various 
cosmological models. 

The linear power spectrum P; is calculated from time integration of the coupled, linearized 
Einstein, Boltzmann, and fluid equations for the metric and density perturbations. The shape 
and growth of the power spectrum depend on cosmological parameters such as the total matter 
density O m , the neutrino fraction Q u , and the Hubble constant h. Simple analytical functions 
were presented that can approximate the linear P; in both CDM and C+HDM models with less 
than 10% error. 



An improved approximation with a higher accuracy has since appeared in Ma (1998). 
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The fully nonlinear power spectrum P n i has to be computed from numerical simulations. 
Simulation output for CDM, LCDM, and C+HDM models was presented to illustrate the sensitive 
dependence of structure formation and evolution on cosmological parameters. We also discussed 
the present understanding of the relation between the linear and nonlinear power spectra, and 
proposed a new approximation which applies to CDM as well as C+HDM models. 

The supercomputing time for this work was generously provided by the National Scalable 
Cluster Project at the University of Pennsylvania and the National Center for Supercomputing 
Applications. 
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k (h Mpc- 1 ) Fig. 1 

Fig. 1. — The linear power spectrum of density fluctuations at z = for the standard CDM model 
(top solid), the C+HDM models with fi„ = 0.1,0.2, and 0.3 (lower solid), and the LCDM model 
with Q rn = 0.3 and SIa = 0.7 (dashed). All are COBE-normalized. 
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Fig. 2. — The smoothed matter distribution at four redshifts, z = 5 (upper left), 3 (upper right), 
1.5 (lower left), and (lower right), from an iV-body simulation of the 0„ = 0.1 C+HDM model. 
Each panel is 100x100 Mpc comoving, and shows projection along one axis of the entire simulation 
box. 
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Fig. 3. — The corresponding power spectrum at z = 3, 1.5, and (from bottom up), for the fi„ = 0.1 
C+HDM model shown in Figure 2. The solid curves show the linear P(k) predicted by the linear 
perturbation theory; the triangles show the nonlinear P(k) computed from the iV-body simulation. 
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Fig. 4. — The smoothed matter distribution at z = 1.5 from iV-body simulations of four 
cosmological models: LCDM with tt m = 0.3 and tl A = 0.7, and C+HDM with VL V = 0.1,0.2, 
and 0.3 (clockwise from upper left). Each panel is 100x100 Mpc comoving, and shows projection 
along one axis of the entire simulation box. The same initial random phases were used in all four 
simulations, so structures appear in similar locations in all panels. 
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Fig. 5. — The matter power spectrum at z = 1.5 for the three C+HDM models shown in Figure 4. 
The solid curves show the linear P(k) predicted by the linear perturbation theory; the triangles 
show the nonlinear P(k) computed from the TV-body simulations. 
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Fig. 6. — The separate CDM and HDM power spectra at z = 3 (lower set) and (upper set) for the 
f2„ = 0.1 C+HDM model. The linear and nonlinear P(k) are shown as solid curves and triangles, 
respectively. For each pair of curves, the higher one is for the CDM and the lower one for the HDM. 
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Fig. 7. — The squares show the ratio of the nonlinear and linear A(k) = 47rk 3 P(k) at z = 
from the VL U = 0.2 (left curve) and 0.1 (right curve) C+HDM simulations. The formula proposed 
by Peacock & Dodds (1996) for the pure CDM models leads to large errors when applied to the 
C+HDM models (dotted curves; left for Q, v = 0.2 and right for 0.1). The solid curves show our 
improved analytical approximation from equation (10), where the right-most one is for the pure 
CDM model. (A modified approximation with a higher accuracy than eq. (10) has since appeared 
in Ma (1998).) 



